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Abstract: We study a complex analogue of a Wilson Loop, defined over a complex 
curve, in non-Abelian holomorphic Chern-Simons theory. We obtain a version of the 
Makeenko-Migdal loop equation describing how the expectation value of these Wilson 
Loops varies as one moves around in a holomorphic family of curves. We use this 
to prove (at the level of the integrand) the duality between the twistor Wilson Loop 
and the all- loop planar S-matrix of A/" = 4 super Yang- Mills by showing that, for a 
particular family of curves corresponding to piecewise null polygons in space-time, 
the loop equation reduce to the all-loop extension of the BCFW recursion relations. 
The scattering amplitude may be interpreted in terms of holomorphic linking of 
the curve in twistor space, while the BCFW relations themselves are revealed as a 
holomorphic analogue of skein relations. 
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1. Introduction 

Scattering processes are usually defined for external states that each have some def- 
inite momentum on the mass-shell. Translational invariance implies that such an 
amplitude is a distribution with support only when 



i=l 

In planar gauge theories the external states are cyclically ordered, so we can solve 
this constraint by introducing 'region' or 'affine' momenta Xi via 



The mass-shell conditions pf = constrain adjacent Xi s to be null separated. In this 
way, the momenta and colour-ordering of the scattering amplitude are encoded in a 
piecewise null polygon, as shown on the left of figure 1. 

In planar A/" = 4 super Yang-Mills, it turns out that scattering amplitudes take 
their simplest form when expressed in (momentum) twistor space. (See e.g. [1, 2, 
3] for introductions to twistor theory and [4, 5] for details specific to momentum 
twistors.) In part, this is because (momentum) twistors make manifest the (dual) 
superconformal symmetry [6] of the tree amplitudes and loop integrand [7], whose 
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Figure 1: A depiction of (a planar projection of) a nodal curve in twistor space corre- 
sponding to a piecewise null polygon in space-time. Via the incidence relations, the vertex 
Xi corresponds to the line Xj = {zi^i,Zi) while the node Zi corresponds to the null ray 
through Xi and Xj+i. 

existence reflects the Yangian symmetry of the planar theory in the amplitude sector. 
(See also [8] for very interesting hints of the role of the Yangian in the loop amplitudes 
themselves.) However, a perhaps deeper reason that these amplitudes belong in 
twistor space is because their twistor data is unconstrained. This is because points 
in space-time correspond to complex lines (linearly embedded Riemann spheres) in 
twistor space, while two such points xi and X2 are null separated if and only if their 
corresponding twistor lines Xi and X2 intersect. Thus, to specify an n-sided piecewise 
null polygon in space-time, we can pick n arbitrary twistors and sequentially join them 
up with (complex^) lines. Adjacent lines intersect by construction, so that as well as 
momentum conservation, the mass-shell condition is automatic. Thus the scattering 
data is encoded in a picture that is usually drawn as on the right of figure 1, where 
the twistors Zi are chosen freely. 

The key observation that motivates this paper is that although figure 1 always 
provides an accurate representation of the twistor curve as an abstract Riemann 
surface, the holomorphic embedding of this Riemann surface in a complex three-fold 
can be far more complicated than this picture indicates. Depending on the choice of 
vertices Zi, i.e. depending on the external kinematics, the embedded twistor curve 
can also appear as depicted in figure 2: that is, it appears as a polygonal knof. 

The new picture hints at an intriguing connection between scattering amplitudes 
and a complexification of knot theory known as holomorphic linking. To understand 
this notion, observe by a simple dimension count that while one can always use 
smooth deformations to untangle real two-surfaces in a real six-manifold, it is not 
always possible to untangle complex curves in a complex three-fold if one is only 

^As usual, we draw real lines in three-space to depict complex lines in the complex three- 
dimensional space. 

^Planarity of the colour-ordered amplitude has nothing to do with the embedding of the curve. 
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Figure 2: Depending on tlie location of the vertices, the twistor lines may form a polygonal 
knot. The plane projection of the nodal curve can thus be far more involved than figure 1 
suggests. 

allowed to vary the curve holomorphically. (This restriction provides a second^ jus- 
tification for the use of the half-dimensional representation in these figures.) That 
holomorphic linking in twistor space should encode interesting information about 
gauge theories in space-time was originally envisaged by Atiyah [9] and Penrose [10]. 
Outside the twistor context, holomorphic linking (and its relation to Abelian holo- 
morphic Chern-Simons theory) has been studied by e.g. Khesin and Rosly [11] and 
by Frenkel and Todorov [12]. 

It is the aim of this paper to elucidate and further investigate this link. Consid- 
erable evidence supporting the existence of such a relation comes from a conjecture 
given by Mason and one of the present authors in [13], building on previous work by 
Alday and Maldacena [14] using the AdS/CFT correspondence and by Drummond, 
Korchemsky and Sokatchev [15], and by Brandhuber, Heslop and Travaglini [16] in 
space-time. The conjecture states, firstly, that the complete classical S-matrix of 
TV = 4 super Yang-Mills is equivalent to the expectation value in the QFT defined 
by the holomorphic Chern-Simons functional [17, 18] 



of a certain operator defined over the complex curve C of figure 2. As we show in 
section 2, this operator is a natural complexification of the standard Wilson Loop 
operator, given by the explicit formula 



^The standard reason is that it gives a better pieture of the structure of the nodes; pictures of 
2-surfaces incorrectly suggest that the component curves are tangent at the nodes. 




(1.3) 




(1.4) 
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in terms of the connection (O,l)-form A, where the precise definition of the meromor- 
phic form co and the meaning of the 'path-ordering' symbol P is explained below. The 
expectation value (WR[C])j^^g is thus a natural complexification of the expectation 
value 



)= j VA e'^cs[^] Tr^jPexp j (1.5) 

of a real Wilson loop in real Chern-Simons theory, that famously computes knot 
invariants [19] such as the HOMFLY polynomial (for gauge group U(A^) and R the 
fundamental representation). 

Extrema of (1.3) correspond to holomorphic bundles on^ CP^'^ that (under mild 
global conditions) in turn yield anti self-dual (or 1/2 BPS) solutions of the M = 
4 super Yang-Mills equations on space-time via the Penrose- Ward correspondence. 
Thus, holomorphic Chern-Simons theory corresponds only to the anti self-dual sector 
of M = 4 SYM. However, it is possible to add a term to (4.1) so that the new action 
corresponds to the complete (perturbative) M = 4 theory [20]. The second part of 
the conjecture of [13] is that, when weighted by the exponential of this new action, 
(W[C]) now yields the complete planar S-matrix of A/" = 4 super Yang-Mills, to all 
orders in the 't Hooft coupling and for arbitrary external helicities. 

These conjectures were supported in [13] by perturbative calculations showing 
that, at least for small k and £, the axial gauge Feynman diagrams of (W[C]) are 
in one-to-one correspondence with MHV diagrams for the i-\oop N'^MHV scattering 
amplitudes. The momentum twistor formulation of the MHV diagram formalism 
had previously been given in [21] and it has since been proved [22, 23] that this 
formalism correctly reproduces the all-loop recursion relation for the integrand of 
planar AT = 4 [7, 24]. 

In this paper, after reviewing the construction of the Wilson Loop operator 
for a complex curve in section 2, we study a complex analogue of the Makeenko- 
Migdal loop equations [25], for both the holomorphic Chern-Simons and A/" = 4 SYM 
theories. At least in principle, these equations give a non-perturbative description of 
the behaviour of (W[C]) as C varies holomorphically. 

Now, the loop equations for the expectation value (1.5) of a fundamental Wilson 
Loop in real U(A^) Chern-Simons theory were shown by Cotta-Ramusino, Guadagnini, 
Martellini and Mintchev [26] to yield the skein relations for the HOMFLY polynomial 
- in other words, recursion relations allowing one to reconstruct the knot invariants. 
Quite remarkably, in section 5 we find that (in the planar limit) the loop equations 
for (W[C]) may be reduced to the BCFW relation [27] that recursively constructs 



The twistor superspace CP ' may be defined as the total space of (^0p3(l))®^ where Cip3(l) 
is the dual of the tautological line bundle on CP'^ and 11 reverses the Grassmann parity of the fibres. 
It is a Calabi-Yau supcrmanifold and we have written for the canonical holomorphic section of 
its Berezinian. 
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the classical S-matrix. The classical BCFW recursion relations thus emerge as a 
holomorphic analogue of the skein relations for a knot polynomial! Replacing the 
pure Chern-Simons theory by the full J\f = 4 theory, the loop equations involve a 
new term and reduce instead to the recursion relation of Arkani-Hamed, Bourjaily, 
Cachazo, Caron-Huot and Trnka [7] for the all-loop integrand. (See also [24] for 
tentative extensions of this relation beyond = 4 SYM.) These results prove the 
conjectures of [13], at least at the level of the integrand. 

Working in space-time, Caron-Huot demonstrated in [28] that the expectation 
value in A/" = 4 SYM of a certain extension of a regular Wilson Loop around the 
null polygon on the left of figure 1, together with extra operator insertions at the 
vertices, also obeys the all-loop BCFW recursion relation. This complete definition 
of this Wilson Loop operator and the vertex insertions was not given in [28], but 
it is expected to be fixed by supersymmetry. Our proof that the twistor Wilson 
Loop (with no extra insertions) obeys the same relation helps confirm that these two 
approaches agree. 

There are at least two ways to view these results. On the one hand, accepting 
the duality between scattering amplitudes and Wilson Loops, one may see the loop 
equations as providing an independent derivation of the all-loop BCFW recursion 
relations. However, we prefer to interpret this work as a field theoretic explanation 
for why the duality holds in the first place: the behaviour of a scattering amplitude 
near a factorisation channel is exactly equal to the behaviour of the twistor Wilson 
Loop near a self-intersection. Away from these singularities, both the scattering 
amplitudes and Wilson Loops are determined by analytic continuation of their data. 
For scattering amplitudes this is just the analyticity of the S-matrix as a function 
of on-shell momenta, encoding the familiar notions of crossing symmetry and space- 
time causality, while for Wilson Loops, analyticity in the twistors Zi expresses the 
fact that (W[C]) generates analytic invariants associated to the holomorphic link. 

Finally, although it is often said that BCFW recursion relations for scattering 
amplitudes have little to do with Lagrangians, path integrals or gauge symmetry, our 
work makes it clear that these standard tools of QFT are indeed intimately connected 
to BCFW recursion from the point of view of Wilson Loops. As always, the Wilson 
Loop is just the trace of the holonomy of a connection around a curve, while the 
loop equations themselves are derived from little more than the Schwinger-Dyson 
equations of the appropriate path integral. We view the connection to Lagrangians 
and path integrals as a virtue, because it suggests that BCFW-type techniques can be 
applied to observables beyond scattering amplitudes. Certainly, the loop equations 
we derive are valid for far more general curves than needed for scattering amplitudes, 
and for more general deformations than BCFW. 
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Figure 3: Going beyond tlie application to scattering amplitudes, we may allow the 
component curves to have arbitrary degrees > 1. 

2. Non-Abelian Holomorphic Wilson Loops 

In this section we review the construction of [13] for a complex Wilson Loop defined 
over a nodal curve C, each of whose components is rational. That is, C = Z{T,) 
where Z is a holomorphic map to CP^I^ and the source curve 



with each component rational and Sj+i fl Sj = {pt} (i being counted mod n). For 
scattering amplitudes, we only need to consider the case that Z restricts to a linear 
map on each irreducible component, so that C generically appears as in figure 2. 
However, it costs little to extend this to the case where each component is mapped 
with arbitrary degree di > 1 (shown in figure 3) and we allow this more general 
possibility. These curves may also be thought of as elliptic curves that have been 
pinched n times around a given homology cycle. They are certainly rather exotic 
from the algebro-geometric point of view^, and it would clearly be interesting consider 
the extension to non-degenerate elliptic curves. 

2.1 Holomorphic frames 

Any C°° bundle E — t- CP'^'^ is inevitably holomorphic on restriction to a arbitrary 
complex curve C, because D^\c = trivially for dimensional reasons. Suppose 
first that C is a single, irreducible rational curve. Provided E\c is holomorphically 
trivial^, we can find a gauge transformation h, varying smoothly over C and such 

^However, it is curious that in the real category, one often considers 'tame' knots - i.e., knots 
that are equivalent to an embedded polygon. 

^By a theorem of Birkhoff and Grothendieck, E\c = ® C(ai) with = if E\c is topolog- 
ically trivial. Triviality as a holomorphic bundle means that, in addition, = for each i and 
(perhaps surprisingly) is in fact generic. 



S = S„ U ■ ■ ■ U S2 u s 
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that 

h-^ o {d + A)\^oh= . (2.2) 

Under this gauge transformation, covariantly holomorphic objects become simply 
holomorphic on C, so /i is said to define a holomorphic frame for E\c- It follows 
from (2.2) that h itself obeys 

{d + A)\^.h = (2.3) 

and so is defined up to a gauge transform h hh', where h' must be globally 
holomorphic on C. 

For two points Zq,z & CP'^'^ that lie on C, we now define 

l]{z,zo) = h{z)h~\zo) (2.4) 

to be the particular solution of (2.3) obeying the initial condition that the gauge 
transformation at 2:0 ^ C* is simply the identity. The holomorphic frame U{z, zo) 
defines a map 

\J{z,Zo):EU,^EU (2.5) 

between the fibres of E at zq and z G CP'^'^, that depends on the choice of curve 
C (not unique if deg(C) > 1). It is the natural complex analogue of the parallel 
propagator U {x, xq) between two points x and xq along a real curve 7. For example, 
it follows immediately from the definition (2.4) that 

U(^,^')U(^',^o) = U(^,^o), (2.6) 

so that U(z, zq) concatenates, just as for the parallel propagator along a real curve. 
In particular, U^z, Zo)'U{zo, z) = U(z, z) = 1r, implying that 

\J{zo,z) = \J{z,zo)-\ (2.7) 

Likewise, if we change our choice of connection (O,l)-form by a smooth gauge trans- 
formation 

A^A,=gdg~'+gAg~\ (2.8) 
then since o (^B + A)\^ o h = {gh)~^ o [p + Ag)\^ o [gk] we find that 

\]{z,z^)-^ g{z)\]{z,z^)g-\zo) (2.9) 

again as for a parallel propagator in the real category. 

In the Abelian case with E a line bundle, setting h = e""^ for some smooth 
function on C, the holomorphic frame equation (2.3) becomes d(j) = A\c- This 
equation always has a solution when C is rational, given by 

(j){z) = / uAA, (2.10) 
Jc 
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where A is restricted to the curve and where w is a meromorphic 1-form on C 
whose only singularities are simple poles at Zq and z, with residues ±1, respectively. 
Exphcitly, if C = Z(S) with Z{a) = z e CP^'^ and Z{ao) = zq, then the pullback of 
a; to E is 

uia')= ^ (2.11) 
^ ^ (a - a'){a' - CTq) 27ri ^ ^ 

in terms of local coordinates a G S. This is just the Green's function for the d- 
operator acting on smooth functions on S. The holomorphic frame equation fixes 
only up to a globally holomorphic function, which by Liouville's theorem must be 
constant. The constant is fixed by the choice of reference point: u vanishes if z and 
zq coincide so (2.10) has (P{zq) = 0. Thus, in the Abelian case we have the explicit 
expression 

U{z,zo) = 6"^ = exp(^- j u AA^ (2.12) 

for our holomorphic frame. 

In the non- Abelian case, acting on a representation R of the gauge group, 11(2, ^o) 
may similarly be given by the (somewhat formal) expression^ 

\J{z,zo) =Pexp j uAA^ . (2.13) 

The meaning of the 'path ordering' symbol P on our complex curve is that the i^^ 
power of the meromorphic differential u that appears from expanding the exponential 
in (2.13) should be taken to have its simple poles located at the reference point zo and 
at either the insertion point of the (i + 1)"^ power of the field A (counting according 
to the colour-ordering), or else at the final evaluation point z. Specifically, 

- c^A^J =1k+ ^(-1)" / /\ (a;(z,) A ^(z,)) (2.14) 

where loi^z^ is a meromorphic differential in Zi with a simple pole at Zq and a simple 
pole at Zj+i (or at z when i = m). Explicitly, in terms of intrinsic coordinates, 

A / N (cr-CTo) dam . dai 

/\uj{ai) = - ^ rir^A---A — . (2.15) 

') [0- - am) ■ ■ ■ {a2 - ai){ai - ao) 27ri 27ri 



''An equivalent formal expression, given in [13], is U(z, zq) = (l + ^^^-^l^) where in the 
series expansion of this expression, the inverse (^-operators are taken to act on everything to their 
right. 
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Equation (2.13) is thus analogous to the (equally formal) expression 



P exp - / A 



X 



XO 



Ir - dy^A,{y) + / dy^A,{y) / dy'^A^iy') - ••• 

J XQ J XQ J XO 

^ f dxf" 

+ 2^ (-!)"" / dsm---dsi Af,{s„^)- — e{sn,-s^-i)---e{s2-si)A^{s 



""-^ ([0,1])" 



dsi 



(2.16) 

for the parallel propagator of a connection d + A = dx'^{dn + A^) along a real curve 
7, where the step functions obey 

j-^e{s - s') = 6{s ~ s') (2.17) 

and are the Green's function for the exterior differential on the interval [0, 1]. We also 
note that the iterated integrals appearing in (2.13) are somewhat similar to iterated 
integrals that arise in studying parallel transport via the Knizhnik-Zamolodchikov 
connection over the configuration space of n-points on a rational curve. See [29] for 
a recent discussion of this in a twistor context. 

2.2 The holonomy around a nodal curve 

Having identified the holomorphic frame U{z,zo) = h{z)h~^{zo) as the natural ana- 
logue of a parallel propagator for a complex curve, we can define the holonomy of the 
connection (O,l)-form around C, based at some point z G CP^'^, to be the ordered 
product 

Hol,[C] = V{Z, Zn) U(^,, Zn-l) " " " U(^2, ^i) U(^i, z) , (2.18) 

where Zi is the image of the node S^+i fl Sj and where, without any essential loss 
of generality, we have taken the base point to be on Z'(Si). Reading from the right 
of (2.18), this definition uses the holomorphic frame \i{zi,z) on ZiT^i) to map the 
fibre -El 2 to the fibre E\z-^ at the node Zi. We can identify this fibre with a fibre of 
E restricted to the next component line Z(T,2), and then map to the following node 
using the holomorphic frame U{z2, Zi) on Z{Tj2)- We continue transporting our fibre 
around the nodal curve in this way until the final factor of U{z, zi) returns us to the 
starting point. The result is thus an automorphism of the fibre E\z, as required. 
The holomorphic Wilson Loop is now defined as in the real category as 

Wn[C] = TTniioUC], (2.19) 

where to take the trace in a representation R we have used the fact that a holomorphic 
frame for E\c induces a holomorphic frame for all tensor products oi E\c and its 
dual, in other words for all choices of representation. Cyclicity of the trace and the 
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concatenation property (2.6) immediately shows that the Wilson Loop is independent 
of the base point. 

For an Abelian gauge group {E a line bundle) the holomorphic frame (2.12) gives 

W[C] =exp(^- J CO A (2.20) 

where w is a meromorphic form on C whose only singularities are simple poles at 
each of the nodes, with residues ±1 on each component. Such a differential is nothing 
but the holomorphic differential ^ on a non-singular elliptic curve in the limit that 
this curve degenerates to our nodal curve (see e.g. [30]). Thus, in the Abelian case 
we recover the prescription 

W[C] =exp(^- j OA (2.21) 

that was proposed by Thomas [31], by Khesin and Rosly [11] and by Frenkel and 
Todorov [12]. For a non- Abelian gauge group we can likewise use (2.13) to formally 
write 

Wr[C] = TiR Pexp(^- j CO A Aj (2.22) 
as a natural complexification of the familiar expression 

^^^[7] = TiR P exp ^ (2.23) 

for a Wilson Loop of a connection on a real curve 7. In equation (2.22), as in (2.13), 
the path ordering symbol means that for each component of C, successive powers of 
u in the series expansion of (2.22) should have one simple pole at one of the nodes 
and another at the subsequent (with respect to the colour-ordering) insertion point 
of the field. 



3. A Holomorphic Family of Curves 

We now investigate the behaviour of our complex Wilson Loop as the curve C varies. 
We shall prove that, just as the change in the real Wilson Loop (2.23) under a smooth 
deformation of a real curve 7 is 

SWr[j] = - dx'' Sx-" TiR (^F^,{x) HoU[7]) , (3-1) 

so too the complex Wilson Loop (2.19) obeys 

SWr[C] = - I uj{z) A dz^ A M Txr{t^-p{z,-z) Ho1,[C]) (3.2) 
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as the complex curve C moves over a holomorphic family. In this equation, 5 is the d- 
operator on the parameter space of the family and, once again, u:{z) is a meromorphic 
differential on C whose only singularities are simple poles at the nodes with residues 
±1. Like its real analogue, equation (3.2) expresses the change in the Wilson Loop 
in terms of the curvature (0,2)-form J-" = dA + A/\A through the complex 2-surface 
swept out by the varying curve. The key point here is that because the curve C varies 
holomorphically, only the (0,2)-form part of the curvature arises. In particular, if E 
is a holomorphic bundle on twistor space, then the Wilson Loop is also holomorphic. 
As before, we prove (3.2) in more generality than we actually need for scattering 
amplitudes, allowing each irreducible component Ci = Z{T,i) to have degree di > 1. 

Here is the proof^. For a fixed partial connection, the holomorphic frame U{zi, zq) 
depends on the choice of two points {zo,zi} G CP^'^ together with a choice of (ir- 
reducible, rational) curve C that joins them. (Note that this curve is not unique 
if its degree is greater than 1.) Equivalently, we can think of this data in terms of 
an abstract Riemann sphere S together with two marked points {aQ,ai} G S and 
a degree d holomorphic map Z : S — CP^'^, such that Z(T,) = C and Z{(Ti) = Zi. 
Now suppose we have a holomorphic family of such maps, parametrized by a mod- 
uli space ^ described, say, by local holomorphic coordinates f. Writing U(cr, (Tq; t) 
for the pullback of U(2;, Zq) to the abstract curve E, for our family of curves the 
holomorphic frame equation is 

{d + Z*A)l]{a, (To; t)=dai — + ^a^(^, t) j U((t, ao; t) = (3.3) 

where we emphasise that both the holomorphic frame and the connection depend on 
t - i.e., on the choice of curve. 

Consider the integral 

uiflia) A U((Ti, a; t) {d + Z*A)V{a, ctq; t) (3.4) 



s 

over S, where Uifi{a) is our familiar meromorphic differential (2.11) whose only 
singularities are simple poles at and (Xi, with residues -|-1 and —1 respectively, 
and where 9 + ^ is as given in (3.3). In principle, the integral depends on these two 
points and on the map Z - that is, it depends on the moduli space - but of course 
it actually vanishes identically on ^ since U((T, (Tq) obeys (3.3). 

*It is also possible to show (3.2) by working (iircctly with the formal expression (2.22), provided 
(as in the real case) due care is taken of the 'path ordering'. 
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Letting 6 be the 9-operator on ^ we have 

= 6 j wi,o(a) A U(ai, a- 1) [d + Z*A) U(a, ctq; t) 
= [ u;i,o(a) AU(cTi,(T;t) (a + ZM)5U(a,ao;t) 



(3.5) 

- / u;i,o(a)U(ai,(7;t) A5(ZM)U((T,(To;t) 



s 



= -5U(ai, ao; t) - ^ u;i,o(a) U(ai, a; t) A 5 (ZM) U(a, ao; t) , 

where the third hne follows from integrating by parts in the first term, then using the 
poles of Uifi and the fact that dU{ai, a; t) = U(cri, a; t)Z*A to perform the integral^. 
In the final term of (3.5) we have 

HZ-A)^\3,A.„ + A.^Ufr.<l,. (S3) 



The most important property here is that, because the map Z varies holomorphically, 
only antiholomorphic derivatives of the connection arise. Again integrating by parts 
we find that 

D<js'U{o-i,ao]t) = 5U(cri,cro;t) + A{z{ai,t))\J{ai,ao;t) - \J{ai,ao;t) A{z{ao,t)) 

ujifl{(j) A da A df U((Ti, a; t) J^^^ — -^U{a, ctq; t) 

uifliz) A d2° A 6z~^ U(2i, z)T^^{z) V{z, zo) 

(3.7) 

where is the natural covariant 9-operator acting on sections of the sheaf ^ over 
^ induced by^° E. 

Equation (3.7) tells us how the holomorphic frame 1]{zi,zq) on an irreducible 
curve varies as we move around in the moduli space. For the Wilson Loop (2.19) on 



c 



In fact, the moduli space may be identified with the Kontsevich moduli space 
of stable, 2-pointed, degree d maps. Because it depends on three points {(To,o'i,<t} e S as well as 
Z, the integrand in (3.4) really lives on the universal curve = Mo,3(CP^''',d) and the pullback 
by the map Z is more properly written as the pullback by the evaluation map on the third marked 
point (labelled a). The integral itself is really the pushdown by the forgetful map -k -.^ ^ and in 
commuting b past the integral sign in (3.5), we really mean the pullback 7r*(5 = d<g — d<fgiss = d^—dn- 
Then (3.6) is just the statement that the evaluation map is holomorphic and so commutes with 5-^, 
with the 9x; term left over. 

""^"More accurately, S = n^.ev'^E. 
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a nodal curve we now readily find 

n 

5W^[C(t)] = Tri? ((:D.^U(z,+i, z.))u(^„ ^.-i) ■ ■ ■ U(z,+2, ^.+i)) 

(3.8) 

= - / ^{z) A dz° A M TrR (j-^^(2, z) Hol,[C(t)]) 

where the integral is taken over the complete nodal curve C{t) = Z{I1) and we 
emphasise that this curve depends on the moduli of the map. As usual, co{z) is a 
meromorphic differential on C{t) with only simple poles at the nodes. This completes 
the proof. 



4. Holomorphic Loop Equations 

In this section we derive an analogue of the loop equations of Makeenko and Migdal [25] 
for our complex Wilson Loops. To do this, instead of considering the complexified 
Wilson Loop itself, we must study its average (Wr[C]) over the space of gauge- 
inequivalent connection (O,l)-forms on E, weighted by the exponential of an appro- 
priate action. 

In section 4.1 we choose the action to be the U(A^) holomorphic Chern-Simons 
functional 

Stcs[A] = ^ J D^I^Atr Q^A^^+^^A^A^^ , (4.1) 

that may be interpreted as the open string field theory of the perturbative B- 
model [17, 18]. In real Chern-Simons theory on a real 3- manifold M, the correlation 
function of a Wilson Loop gives a knot invariant. Since these invariants depend 
only on the (regular) isotopy class of the knot, one expects that, unlike W/j[7] itself, 
(W/{[7])(-,g should remain invariant as one smoothly deforms 7. This is essentially 
true, except that there are important corrections when the deformation causes 7 to 
pass through a self-intersection so that the isotopy class of the knot jumps. The 
loop equations determine the behaviour of (WR[7])^g under such a jump and it was 
shown in [26] that, for Wilson Loops in the fundamental representation of U(A^), 
they amount to the skein relations 

s^P+(r, s) — s^^P^{r, s) = (r^ — r^^)Po{r, s) (4.2) 

from which the HOMFLY polynomial P-y{r, s) can be recursively constructed^^. 

^^The variables of the HOMFLY polynomial are related to N and the level k of the Chern-Simons 
theory by and s = and r = e"^/^, where A = 2'KiN/{k + N) may be interpreted as the 't Hooft 
coupling. In [26], the skein relations were obtained to lowest order in A. The skein relation only 
determines the HOMFLY polynomial for Wilson Loops in the fundamental. 
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We expect that correlation functions of our Wilson Loops in the holomorphic 
Chern-Simons QFT are likewise associated with holomorphic linking invariants. In- 
deed, [11, 12] show that for two non-intersecting curves Ci, C2 C of genera gi > 1, 
the expectation value (W[Ci] W[C2]) in Abelian holomorphic Chern-Simons theory 
on involves^^ the holomorphic linking invariant 

L{Cue,;C2;e2) = — / ^ d^^" A dw;^ A A ^2 (4.3) 

Ait J ■' \z — 

C1XC2 

that generalizes the Gauss linking number 

i^(7.,7,) = ^ / e,,,j^^dx^ Ad/ (4.4) 
Air J \x — y]-^ 

in for real Abelian Wilson Loops in M.^. Observe that (4.3) depends on a choice of 
holomorphic 1-form 6i on each curve. For our nodal curve, the corresponding form is 
meromorphic and depends on the location of the nodes. In section 4.1 we shall find 
that (5(Wi^[C])j^|^g generically vanishes, with important corrections when the type 
of the curve jumps. In the planar limit, the loop equations may be reduced to the 
BCFW relations [27] that recursively construct the classical S-matrix. In this sense, 
the BCFW relations are simply the U(oo) skein relations for the holomorphic link 
defined by C. 

As mentioned in the introduction, holomorphic Chern-Simons theory on CP^'^ 
corresponds only to the anti self-dual sector of A/" = 4 super Yang-Mills theory on 
space-time. In section 4.2 we replace (4.1) by the action [20] 

Sm=4A] = 5hcs[^] + j d^l'x lndet(a + ^)|3, (4.5) 

that, as our notation indicates, describes A/" = 4 super Yang-Mills theory on twistor 
space. We shall see that the new term is itself eminently compatible with the Wilson 
Loop. The loop equations in this theory involve a new term that leads to the all-loop 
extension of the BCFW relations found by Arkani-Hamed et al. in [7]. 

4.1 Holomorphic linking 

In this section, we choose the action to be the holomorphic Chern-Simons func- 
tional (4.1), where D^^^z is the canonical holomorphic section of the Berezinian of 
CP^'^, given explicitly by 

D=^IS = ^eai3ysz"dz'^ A dz^ A dz^ ^eabcddip" dip'' dip" dip'^ (4.6) 

-'^^We will sec later that, at least for the Af — 4: theory, 'holoniorphic self-linking' is well-defined 
even without a choice of framing. 
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in terms of homogeneous coordinates (z°',ip"') E C'*'^. This action is defined on the 
space of partial connections Z) on a C°° vector bundle E — )■ CP^'^, given at least 
locally hj D = d + Ain terms of a background partial connection d. We take this to 
obey 3^ = and so define a background complex structure on E. Thus, A{z, z, ip) 
is a superfield with component expansion 

A{z, z, ^) = a{z, z) + V"r,(z) + ^^>^ ^ah{z, z) + --- + ^eabcdri^'ri^" 9{z, z) , 

(4.7) 

'■^ valued in smooth sections of 



where the coefficient of {ipY is a (O,l)-form on 

End(E) ® C>cp3(-^)- 

We choose the gauge group to be U(A^) and consider only Wilson Loops in the 
fundamental representation. The Yang-Mills (or open string) coupling g^ is then 
related to the Chern-Simons level by g^ = 27r/(A; + A^). In anticipation of taking 
the planar limit, we write = \/N in terms of the 't Hooft coupling A and include 
a factor of in definition of the Wilson Loop so that W[C] = 1 if the holonomy 
is trivial. 



Inserting (3.2) into the holomorphic Chern-Simons path integral, one finds 
1 



5(W[C(t)]> 



N 
X 



VA 
VA 



U[Z] 



C{t) 



A d^° A Szl^ tr( J^,^(z) HoL[C(t)] 



■Shcs[^] 



C(t) 



uj{z) At! Hol,[C(t)] 



6 



5A{z) 



~Shcs[^] 



(4.8) 

since the variation of the holomorphic Chern-Simons functional is J^'"^ x N/X. To 
obtain an interesting equation for 5(W[C(t)]), as in [25] we integrate by parts in 
the path integral, bringing the variation of the connection to act on the holonomy. 
Essentially the same argument as in section 3 shows that under a variation of the 
connection at some point z' G C, (displaying colour indices) 

ui,o{z') A U(zi, z') {d + A)^\J{z', zo) ^ 







6 



5Aiz) ^ 
6 



c 



U(zi,zo)^ 



(4.9) 



co^,o{z')A\J{z^,z')'J'^\z,z')\J{z',Zo) 



^Az) , J JC 

where we have used the fact that, for a U(A^) gauge group, 

{6l6A{z))\A{zTn = '^'^\z,z')5\6-^ 

with 5^1^(2,2;') a Dirac current concentrated on the diagonal in CP^'^ 
5^\^{z, z') is a distribution-valued (0,3)-form on CP^I^ X 
fiO'PrCP^I^) we have 



(4.10) 

,3|4 . 



i3|4 



CT3|4 



) such that for any a G 
T)''\^z' A 6^\\z, z') A a{z') = a{z) . (4.11) 
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Figure 4: The expectation value ^W[C(t)]) in holomorphic Chern-Simons theory varies 
holomorphically except where the curve develops a new node. 



Explicitly, S^^'^{z,z') may be represented by the integral [5, 32] 

= / ^ A mz - uz') = / ^ A ^i?^ tii'i'' - -*'°) (4-12) 

a=l a=l 

that forces oc z'"' so that the two point must coincide projectively. 

We now apply equation (4.9) to the variation holonomy based, say, at a point 
z E Ci. One finds 

tr (^^HoL[C(t)]) = [u(^,^„)\U(^.,^._i)V--U(^i,^r; 

= V / C0i+^,^{z') A z') tr \\J{z, z^) ■ ■ ■ \J{zi+i, z')] tr \\J{z', z^) ■ ■ ■ U(zi, z) . 

i Jc^t) L J L J 

(4.13) 

The factor of S^^^{z,z') ensures that this term has support only at points t in the 
moduli space where the curve C{t) degenerates so that the component Ci(t) ei- 
ther self-intersects or else intersects some other component. Similar contributions 
arise from summing over the possible locations of the base-point z. In other words, 
thinking of our original curve as an elliptic curve, (5(W[C(t)]) vanishes everywhere 
except on boundary components of the moduli space where the holomorphic map 
degenerates so that two points on the source curve are mapped to the same image. 

Combining all the terms, the holomorphic Chern-Simons loop equations (4.8) 
become 

6{W[C{t)]) = -X J uj{z) Auj{z') Ap\\z,z') {W[C'{t)]W[C"{t)]) , (4.14) 

C(t)xCit) 

for normalised Wilson Loops in the fundamental of a U(A^) gauge group, where C'(t) 
and C"{t) are the two curves obtained by ungluing C{t) at its new node z = z' (see 
figure 4). In the planar limit, the correlation function of a product of Wilson Loops 
factorizes into the product of correlation functions, so (4.14) simplifies to 

S{W[C{t)]) = j ^(^) A uj{z') A 5'\\z, z') (W[C"(t)]) (W[C"(t)]) (4.15) 

C{t)xC[t) 
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where A is the 't Hooft couphng. Thus, in the planar hmit, the behaviour at a self- 
intersection is determined by the product of the expectation values of the Wilson 
Loops around the two lower degree curves. This may be viewed as an analogue 
for holomorphic linking of the skein relations for real knot invariants. In section 5 
we shall see that the classical BCFW recursion relations are simply a special case 
of (4.15). 

4.2 Loop equations for the twistor representation oi J\f = A SYM 

We now consider the expectation value of the complex Wilson Loop in the twistor 
action 

Sm=4 [A] = ^hcs [A] + d^l^x log det {d + A)\^ , (4. 16) 

where X is a line in twistor space and F is a middle-dimensional contour in the 
space Gr2(C^) of these lines^'^, interpreted as a choice of real slice of complexified, 
conformally compactified space-time. The proof that (4.16) is an action for A/" = 4 
SYM on twistor space is given in [20], to which the reader is referred for further 
details. (Readers familiar with the MHV formalism may note that the new term 
in (4.16) provides an infinite series of MHV vertices, while the 3-point MHV vertex 
in the Chern- Simons theory vanishes in an axial gauge. From the point of view 
of twistor string theory [18], the determinant may be understood as the partition 
function of a chiral free fermion CFT on X, coupled to the puUback of the twistor 
space gauge field.) 

From the point of view of the present article, it is possible to motivate this action 
by observing that Wilson Loops on complex curves are again the natural observables 
of the QFT based on (4.16). The reason this is so is because of the close connection 
between determinants and holonomies. More precisely, just as 

5 log det M = 5trlogM = ii{M~^5M) (4.17) 

for a finite dimensional matrix M, Quillen shows [33] that, for an abstract curve S, 
under a variation of the connection, the logarithm of a section of the determinant 
line bundle Det — )■ s^^'^ varies as^^ 

5\ogdei{d + A)\^= [ tT {J_A A 6A) , (4.18) 

where 

JA{a) = lim (CAia, a) - Go{a' , a)) (4.19) 
^^The fermionic integral is done algebraically, as always. 

^^In writing this formula, we assume that 5 is a flat connection on Dot in a trivialisation where 
the local connection 1-form vanishes. That Det should admit such a flat connection is necessary if 
we wish to treat det(9 + A)y: as a (holomorphic) function on (the space of connection (0,1)- 

forms on S) . However [33] , the construction of this flat connection requires a choice of background 
connection. 
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is the limit^^ on the diagonal in S x S of difference between the Green's function 
for ^ + and the Green's function for background connection d^, ■ These Green's 
functions are non-local operators on the Riemann surface; the fact that we take their 
limit on the diagonal can be understood as part of the trace. For a single such Green's 
function, the limit on the diagonal is necessarily singular, but the singularity cancels 
in the difference (4.19). The dependence on a background connection reflects the fact 
that to treat sections of the determinant line bundle as functions we must first pick 
a trivialisation, amounting to a choice of background connection (O,l)-form on the 
bundle over S. In our case, we will of course choose this background (9-operator to 
be the one induced by our choice of base-point d G ^■^'^ in writing the holomorphic 
Chern-Simons action. 

When E is a Riemann sphere, the Green's functions are 

Gj,{a', a) = h{a')Go{a', a)h'\a) and Go(a', a) = ' (4-20) 

ZTTl a — (J 

in terms of the holomorphic frame h on E. Therefore Quillen's prescription reduces 
to 

da U(ct',ct) -U(ct,(t) da d\]{a',a) 



Ja{(^) = — r lim 



27ri a'^a a' — a 27ri da' 



(4.21) 



For our purposes, the presence of the holomorphic derivative here is somewhat in- 
convenient. However, it is easily seen from the formal series (2.14) that l]{a',a) 
depends only meromorphically on a' and is regular at a' = a. Thus we can use 
Cauchy's theorem to rewrite the derivative in (4.21) as an integral, so 

5 log det {d + A)\^ = jdaA tr (^^ ^-^^U(a', a) 5A{a)^ 

s 

1 f dcr A da' A do-" , 

-tr(U(a',a)5^(a)) 



(27ri)3 J (a - a') {a' - a") {a" - a) 

(4.22) 

where the integrals over a" and a' are performed over contours encircling the poles at 
a" = a' and a' = a, respectively. (The reason for introducing the new point a" G E 
will become clear momentarily.) 

To apply this to the case that E = Z~^(X) for X a line in twistor space, suppose 

Z{a") = Za Z{a') = Zb and Z{a) = z = Za + az^ (4.23) 

so that X = Span[2;a, Zb\. Then 

da A da' A da" 



(a - a'){a' - a"){a" - a) 



■ (ada) A(&d6) 

A 



(4.24) 



""^^Below, this limit will be revealed as nothing but the forward limit of an amplitude! 
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where u}a,b{z) is meromorphic with only simple poles only z = Za and z^, and where 
{ada){b db) / (ab)^ is the fundamental bi-differential on X with a double pole along 
the diagonal. This bi-differential combines with the integral over the space of lines 
X, since 

d^l^x A = D^|4,^ ^ , (4.25) 

(a b)^ 

giving a contour integral over CP^I^ X CPf . Therefore, the variation of the new term 
in the action is 



S J d^l^a; logdet(9 + ^) 



(27ri 

r rxxxsixsi 



X- / (0^\\3 AuJa,biz)tl{\]{Zb,z) 5A{z) 

(4.26) 

where the 5*^ x contour sets Za ^ z and Zb — )■ z, differentiating \J{zh, z) in the 
process. The integral over all i G X is performed using the dependence on z in the 
holomorphic frame and in the variation of the connection. 

Equation (4.26) expresses the variation of the logarithm of the determinant in 
the twistor action (4.16) in terms of a holomorphic frame on the line X. It dovetails 
beautifully with the holonomy around the curve C, providing a new contribution 



A 



u{z)AUa,b{z)AS''\\z,z) tT(V{Zb,z) Hol,[C(t)] 

rx5ix5i \Cit)xX 

(4.27) 

to the loop equations. The 5-function 5^^^{z,z) in this expression means that this 
term only contributes at points in the moduli space where the curve C{t) 3 z inter- 
sects the line X 3 z. The integral over F then adds up these contributions for each 

X G r. _ 

To interpret the product 11(2;;,, z) B.o\z[C(t)], replace C{t) by a curve C{t) where 

C(t)nX= {z, Zb} G CP=^l^ (4.28) 

and such that as Zb — )■ z, C(t) C(t). As a nodal curve, C{t) has one more 
component than C{t), with one component becoming double covered in the limit 
(see figure 5). Now, because the contour integral forces Zb = z, we can replace the 
holonomy around C(t), based at 2; = f = C(t) fl X, by a product of holomorphic 
frames that transport us around C{t) from Zb to z. The remaining factor of U{zb, z) 
transports this frame back to Zb along X. The resulting trace can be interpreted as 
Wilson Loop around C {t) U X, so that 

^tr(u(zfe,z)Hol,[C(t)]) = W[C^) UX] , (4.29) 
provided we are inside the contour integral over the location of Zb C X. 



-19- 



X 



c 



z = z 



X 



c 



Figure 5: The new term in the loop equations for = 4 SYM can be understood as a 
holomorphic Wilson Loop around a new curve C U X, where C reduces to C when — )• z 
along X. Note that the nodal curve C has one more component than C . In the scattering 
amplitude context, (W[C U X]) may be interpreted as a forward limit of an n + 2 particle 
amplitude. 



Finally then, we see that the holomorphic loop equations for the twistor repre- 
sentation of planar A/" = 4 SYM are 

5(W[C(t)]) = - A y u{z) A A b''\\z, z') (W[C"(t)]) (W[C"(t)]) 

C(t)xC(t) 



A y D^IS, A D=^IS, 
rxS'lxS'i 



uj{z) A UaA^) A z) (W[C(t) U X]> 



C(i)xX 



(4.30) 

and, as in the original equations of Migdal & Makeenko [25], are entirely expressed 
in terms of expectation values of our (complexified) Wilson Loops. In the following 
section we show that for a particular one-parameter family of curves C(t), this holo- 
morphic loop equation reduces to the all-loop BCFW recursion relation found in [7], 
but we emphasise that (4.30) is valid for a far more general class of curves. 



5. BCFW Recursion from the Loop Equations 

We now apply the holomorphic loop equation (4.30) to the case that each component 
of C is a line, as in figure 2. C can be completely specified by giving the location 
of its n nodes, and we shall often write C = (12) U (23) U ■ ■ ■ U (nl) for the curve 
(where = Span[2;j, Zi+i]) and W[l, 2, . . . ,n] for the corresponding Wilson Loop 

(in the fundamental representation). We wish to show that the loop equations for 
(W[l, 2,...,n]) directly reduce to the all-loop BCFW recursion relations of Arkani- 
Hamed et al. [7], although it is important that our derivation of these relation is 
completely independent of any notions of scattering theory. 
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Figure 6: The family of nodal curves C(t) corresponding to a BCFW deformation of a 
planar amplitude in A/" = 4 SYM, shown for n = 8. As we vary over the moduli space, the 
line (l,n(t)) sweeps out a pencil of lines in the plane (n— l,n, 1) with focus zi. This plane 
necessarily intersects every other component of C{t). 

Consider the holomorphic family of Wilson Loops determined by translating the 
node Zn along the line (n — 1,?t,). That is, we define 

Zn(t) = Zn + tZn-l (5-1) 

and consider the 1-parameter family of holomorphic curves 

C{t) = (12) U (23) U ■■■ U (n-ln(t)) U (n(t) 1) (5.2) 

as shown in figure 6. Notice that although t can be thought of as a local holomorphic 
coordinate on the line (n— Ira), it is really a coordinate on the moduli space of our 
family of curves, i.e. a coordinate on a CP^ C As t varies, the line (ra(t), 1) sweeps 
out the plane (n — l,n,l) with all other components of C(t) remaining fixed. Lines 
and planes necessarily intersect 16 in CP^ so for every j = 3, . . . , n—1 there exists a tj 
for which {n{tj), 1) intersects the component Cj. Let us call this intersection point 
Ij and write Uj for the point Zn(tj). Clearly, we have 

= (n-l,n, 1) n (j-1, j) and % = (n-1, n) fl (j-1, j, 1) , (5.3) 

essentially by definition. (Note that we do not get any new intersections when j = 
n,l,2.) 



^^Lines and planes do not necessarily intersect in the superspace CP^'"', because they might 'miss' 
in the fcrmionic directions. The fermionic 5-functions in the S^^'^{z,z') ensure that the right hand 
side has support only when the curves do intersect in the superspace. See also the discussion in [21]. 
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We now study the loop equations for this family. Considering first the case of 
pure holomorphic Chern-Simons theory, (4.15) becomes 



' d/ 

-A6{W[C{t)]) = 

A J2 I jAujiz)Au;{z')Ap\\z,z') {W[C'{t)]) {W[C"{t)]) 



CF'^xCi(t)xCj 



where 

C"(t) = (12)U(23)U---U(j-l,/,) 
C"(t) = (J,,j)U(j,j + l)U---(n,-,l) 

and where we have integrated (4.15) over our moduli space CP^ C ^ using the 
meromorphic differential dt/t. We can dispose of the left hand side immediately: 
integrating by parts gives 

^ A 6{W[C{t)]) = {W[Cm) - (W[C(oo)]> 

I (5.0) 
= (W[l,2,...,n]>-(W[l,2,...,n-l]>, 

the difference of the Wilson Loops around the original curve and the curve with 
Zn = Zn-i- Notice that the holomorphic frame U(i„, Zn-i) — t- 1 as — ^ -^n-i so that 
while the line {n—l,n) does not simply 'disappear' from the picture, neither does it 
contribute to the holonomy. 

The right hand side of (5.4) is almost as straightforward to compute, because the 
integrals over both the CP^ moduli space and the two copies of C (t) are completely 
fixed by the 5-function S^^^(z,z'). Since z G (n(t),l) and z' G (j — l,j) we can 
parametrize these integrals by setting 

z = Zn(t) + szi = Zn + + szi and z' = Zj_i + rzj (5.7) 

in terms of local coordinates r and s on Ci{t) and Cj, respectively. With these 
coodinates, the meromorphic forms uj{z) and oo{z') become simply 

d dv 

u{z) = — and u{z') = — (5.8) 
s r 

which indeed have simple poles at the nodes of Ci (t) and Cj as required. Performing 
the integral is then merely a matter of using the explicit form (4.12) of 5^^^{z, z') and 
computing a Jacobian. In fact, this integral is 

r 1 [ drdsdt du -^.u, . ,_ ^, 

n-l,ra,l,j-l,j = / b ^ [Zn + tZn-i + szi+uzj^i+rzj) 5.9) 

J r s t u 
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which was shown in [5] to be just the basic dual superconformal invariant Ri-jn of [6]. 
Including the product of the two Wilson Loops on the smaller curves we find that 
the loop equations in pure holomorphic Chern-Simons theory reduce to 

(W[l,...,n]) = (W[l,...,n-l]) 

n-1 

+ A ^ [n-l,n,l,j-l,j](W[l,...,j-l,/,])(W[/„j, ...,n-l,%]), 

i=3 

(5.10) 

where Ij and nj were given in (5.3). Re-interpreting the twistor space as momentum 
twistor space, this is just the tree-level BCFW recursion relation [27] in the momen- 
tum twistor form given in [7], summed over all MHV degrees. Expanding (5.10) in 
powers of the Grassmann coordinates at each vertex gives the BCFW relations for 
specific N'^MHV partial amplitudes. 

Turning now to the full JV = 4 theory, from (4.30) we have an additional contri- 
bution 

rxSix5i [cpixC(i)xX 

(5.11) 

where 

C(t)UX = (12)U---(n-l,n(t))U(z,6)U(6, 1) (5.12) 

is our n + 2 component curve. Once again, the integrals inside the square brackets 
are completely fixed by the P^^{z, z). The intersection point z and the point h{t) are 
fixed to be 

^ = (a, 6) n (n — 1, n, 1) and n(t) = (n— 1, n) fl (a, 6, 1) (5.13) 

and the Jacobian from the t, z and z integrals gives the /^-invariant [n — 1, ra, 1, a, h]. 
The N' = 4 loop equations thus reduce to 

(W[l,...,n]> = (W[l,...,n-l]> 

n-1 

+ A 5^ [n-l,n,l,j-l,j](W[l,...,j-l,/,]>(W[/„j,...,n-l,n,]) 

i=3 

+ A y D^lS.AD^lSfe [n-l,n,l,a,6](W[l,...,n-l,n,b,l,26]>, 

rx^ixsi 

(5.14) 

where we remind the reader that the x 5*^ contour is taken to fix Za^ — ?■ z along X 
(the i?-invariant has a simple zero in this limit, cancelling one of the factors of {za z^) 
in the denominator of the measure). Re-interpreting the twistor space as momentum 
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twistor space, this is the extension of the BCFW recursion for the all- loop integrand 
found by Arkani-Hamed et al. in [7]. 

Let us conclude with a few remarks. Firstly, the explicit power of the 't Hooft 
coupling A is present in (5.14) because our normalisation of the twistor action cor- 
responds to the normalisation ^ J tr(F A *F) = ^ / tr(F A *F) of the space-time 
Yang-Mills action. In this normalisation, n-particle, £-loop scattering amplitudes are 
proportional to X^~^N^ (with x = 2 for planar diagrams). If one rescales the twistor 
(or space-time) connection as A ^ \/XA, the explicit A disappears from (5.14) and 
the £-loop planar scattering amplitude corresponds to the coefficient of ^("■+2^-2)/2 

Secondly, when computing real knot (rather than link) invariants, it is necessary 
to choose a 'framing' - roughly, a thickening of the knot into a band. For Wilson 
Loops in real Chern-Simons theory, this may be understood perturbatively as a 
point-splitting procedure that regularises potential divergences as the two ends of a 
propagator approach each other along the knot 7. This would seem to be particularly 
important in the case that 7 has cusps. Above, however, we made no attempt to 
choose a 'holomorphic framing' of C. In fact, at least for the pure holomorphic 
Chern-Simons theory, this does not seem to be necessary. Consider for example a 
propagator stretched between two adjacent line components of C, say — and 

+ We can regularise any potential divergence in this propagator by ungluing 
the node, replacing + — + where = Zi + ez' for some arbitrary 
point z'. As shown in [13], the integral of an axial gauge holomorphic Chern-Simons 
propagator along two non-intersecting lines gives an i?-invariant, and in this case we 
obtain — + (where 2* determines the axial gauge). It is then easy to 

show that this is of order e as e — )■ 0, as a consequence of the A/" = 4 supersymmetry. 
In other words, the short- distance singularity of the propagator is always integrable, 
provided one sums over the complete A/" = 4 supermultiplet. 

If the holomorphic Chern-Simons expectation value is finite (for generic Zi), 
the same cannot be said of the expectation value in the full theory. The integral 
over the space of lines X will diverges (at least if F is a contour corresponding 
to a real slice of space-time), with the divergence coming from those X G F that 
intersect C at more than one place, or if the intersection occurs at a node, and 
refiecting the divergence of the corresponding scattering amplitudes in the infra- 
red. Consequently, without a regularisation (or else a choice of 'leading singularity' 
contour for F), the expression (5.14) is somewhat formal, and provides a recursion 
relation for the all-orders integrand, rather than the all-orders Wilson Loop itself. It 
would be fascinating if it is possible to regularise these divergences using some form 
of holomorphic framing of C{t). Curiously, the Coulomb branch regulator of [34] - 
interpreted in dual conformal space-time as bringing the vertices of the piece-wise 
null polygon into the interior of AdSs - corresponds in twistor space to replacing 
the line X by a (complex) quadric. See e.g. [35, 36, 37, 38] for discussions of this 
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regulator in a twistor context. 

Further, although we considered only the particular deformation z„ — )■ 
relevant for BCFW recursion, it should be clear that the loop equations can equally 
be applied to the more general deformation [39] 

— )■ — tCjZ* i = l,...,n, (5.15) 

of the scattering amplitude curve, that translates each vertex towards some arbitrary 
point G CP^'^. (Here, q are arbitrary complex constants.) This one-parameter 
deformation includes both BCFW recursion and the MHV diagram formalism as 
special cases and the validity of the corresponding recursion relation for the all-loop 
A/" = 4 integrand has been shown recently in [22]. It is also clear that many other 
deformations of C are possible and it would be interesting to know if anything useful 
can be learned from, say, non-linear or multi-parameter deformations. 

As a final remark, many possible extensions to this work suggested by appealing 
either to the analogy with real Chern-Simons and knot theory, or else to the link 
to scattering amplitudes in space-time. We hope that this close relationship helps 
to fertilise new progress, both in understanding holomorphic linking and in A/" = 4 
SYM. 

Acknowledgments 

We thank Freddy Cachazo, Jaume Gomis and Natalia Saulina for helpful comments. 
We would also particularly like to thank Lionel Mason for many useful discussions. 
The work of DS is supported by the Perimeter Institute for Theoretical Physics. Re- 
search at the Perimeter Institute is supported by the Government of Canada through 
Industry Canada and by the Province of Ontario through the Ministry of Research 
& Innovation. The work MB is supported by an STFC Postgraduate Studentship. 

References 

[1] R. Penrose and W. Rindler, Spinors and Space-Time, vol. 2. Cambridge University 
Press, 1986. 

[2] S. Huggett and P. Tod, An Introduction to Twistor Theory. Student Texts 4. London 
Mathematical Society, 1985. 

[3] R. Ward and R. Wells, Twistor Geometry and Field Theory. CUP, 1990. 

[4] A. Hodges, Eliminating Spurious Poles from Gauge-Theoretic Amplitudes, 
arXiv: 0905. 1473 [hep-th] . 

[5] L. Mason and D. Skinner, Dual Superconformal Invariance, Momentum Twistors 
and Grassmannians, JHEP 11 (2009) 045, [arXiv : 0909 . 0250]. 



-25- 



[6] J. M. Drummond, J. Henn, G. P. Korchemsky, and E. Sokatchev, Dual 

Superconformal Symmetry of Scattering Amplitudes m = 4 Super-Yang-Mills 
Theory, Nucl. Phys. B828 (2010) 317-374, [arXiv : 0807 . 1095]. 

[7] N. Arkani-Hamed, J. Bourjaily, F. Cachazo, S. Caron-Huot, and J. Trnka, The 
All-Loop Integrand for Scattering Amplitudes in Planar A/" = 4 SYM, 
arXiv: 1008.2958 [hep-th] . 

[8] J. Drummond, L. Ferro, and E. Ragoucy, Yangian Symmetry of Light-Like Wilson 
Loops, arXiv: 1011 .4264. 

[9] M. F. Atiyah, Green's Functions for Self-Dual Four Manifolds, Adv. Math. Supp. 7 A 
(1981) 129-158. 

[10] R. Penrose, Topological QFT and Twistors: Holomorphic Linking, Twistor 
Newsletter 27 (1988) 1-4, [http://people.maths.ox.ac.uk/lmason/Tn/]. 

[11] B. Khesin and A. Rosly, Polar Homology and Holomorphic Bundles, 
Phil. Trans. Roy. Soc. Lond. A359 (2001) 1413-1428, [math/0102152]. 

[12] I. Frenkel and A. Todorov, Complex Counterpart of Chern- Simons- Witten Theory 
and Holomorphic Linking, Adv. Theor. Math. Phys. 11 (2007), no. 4 531-590, 
[math/0502169]. 

[13] L. Mason and D. Skinner, The Complete Planar S-matrix of N = 4 Super Yang-Mills 
from a Wilson Loop in Twistor Space, JHEP 12 (2010) 018, [arXiv : 1009 . 2225]. 

[14] L. F. Alday and J. Maldacena, Gluon Scattering Amplitudes at Strong Coupling, 
JHEP 06 (2007) 064, [arxiv 0705.0303 [hep-th]]. 

[15] J. Drummond, G. Korchemsky, and E. Sokatchev, Conformal Properties of 
Four-Cluon Planar Amplitudes and Wilson Loops, Nucl. Phys. B795 (2007) 
385-408, [arXiv: 0707. 0243 [hep-th]]. 

[16] A. Brandhuber, P. Heslop, and G. Travaglini, MHV Amplitudes in N = A Super 
Yang-Mills and Wilson Loops, Nucl. Phys. B794 (2008) 231-243, 
[arXiv: 0707. 1153]. 

[17] E. Witten, Chern-Simons Cauge Theory as a String Theory, Prog. Math. 133 (1995) 
637-678, [hep-th/9207094]. 

[18] E. Witten, Perturbative Cauge Theory as a String Theory in Twistor Space, 
Commun. Math. Phys. 252 (2004) 189-258, [hep-th/0312171]. 

[19] E. Witten, Quantum Field Theory and the Jones Polynomial, Commun. Math. Phys. 
121 (1989) 351. 

[20] R. Boels, L. Mason, and D. Skinner, Supersymmetric Cauge Theories in Twistor 
Space, JHEP 02 (2007) 014, [hep-th/0604040]. 



-26- 



[21] M. Bullimore, L. Mason, and D. Skinner, MHV Diagrams in Momentum Twistor 
Space, JHEP 12 (2010) 032, [arXiv: 1009. 1854]. 

[22] M. Bullimore, MHV Diagrams from an All-Line Recursion Relation, 
arXiv: 1010.5921. 

[23] S. He and T. McLoughlin, On All-Loop Integrands of Scattering Amplitudes in 
Planar J\f = 4 SYM, arXiv: 1010.6256. 

[24] R. Boels, On BCFW Shifts of Integrands and Integrals, arXiv: 1008.3101 [hep-th] . 

[25] Y. Makeenko and A. Migdal, Exact Equation for the Loop Average in Multicolor 
QCD, Phys. Lett. B88 (1979) 135-137. 

[26] P. Cotta-Ramusino, E. Guadagnini, M. Martellini, and M. Mintchev, Quantum Field 
Theory and Link Invariants, Nucl. Phys. B330 (1990) 557. 

[27] R. Britto, F. Cachazo, B. Feng, and E. Witten, Direct Proof Of Tree-Level Recursion 
Relation In Yang-Mills Theory, Phys. Rev. Lett. 94 (2005) 181602, 
[hep-th/0501052]. 

[28] S. Caron-Huot, Notes on the Scattering Amplitude / Wilson Loop Duality, 
arXiv: 1010. 1167. 

[29] Y. Abe, Holonomies of Gauge Fields in Twistor Space 1: Bialgehra, Supersymmetry 
and Gluon Amplitudes, Nucl. Phys. B825 (2010) 242-267, [arXiv: 0906. 2524]. 

[30] J. Fay, Theta Functions on Riemann Surfaces. Lecture Notes in Mathematics. 
Springer, 1973. 

[31] R. P. Thomas, A Holomorphic Casson Invariant for Calahi-Yau 3-folds, and Bundles 
on K3 Fihrations, J. Diff. Geom. 54 (2000), no. 2 367-438, [math/9806111]. 

[32] L. Mason and D. Skinner, Scattering Amplitudes and BGFW Recursion in Twistor 
Space, JHEP 01 (2010) 064, [arXiv: 0903. 2083]. 

[33] D. Quihen, Determinants of Cauchy- Riemann Operators over Riemann Surfaces, 
Func. Anal. Appl. 19 (1985) 37-41. 

[34] L. F. Alday, J. Henn, J. Plefka, and T. Schuster, Scattering Into the Fifth Dimension 
ofM = 4: Super Yang-Mills, JHEP 01 (2010) 077, [arXiv : 0908 . 0684 hep-th]. 

[35] A. Hodges, The Box Integrals in Momentum Twistor Geometry, arXiv : 1004 . 3323. 

[36] L. Mason and D. Skinner, Amplitudes at Weak Coupling as Polytopes in AdS^, 
Submitted to JPhys A (2010) [arXiv: 1004.3498]. 

[37] L. F. Alday, Some Analytic Results for Two-Loop Scattering Amplitudes, 
arXiv: 1009. 1110. 



-27- 



[38] J. Drummond, J. Henn, and J. Trnka, New Differential Equations for On-Shell Loop 
Integrals, arXiv: 1010 . 3679. 

[39] H. Elvang, D. Z. Freedman, and M. Kiermaier, Proof of the MHV Vertex Expansion 
for All Tree Amplitudes in M = A SYM Theory, JHEP 06 (2009) 068, 
[arXiv: 0811. 3624]. 



-28- 



